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The experimental observable properties of the triplet p z -wave pairing state, proposed by Wu et al. 
[arXiv: 1503.06707] in quasi-one dimensional A 2 Q- 3 AS 3 materials, are theoretically investigated. This pair¬ 
ing state is characterized by the line nodes on the k z = 0 plane on the Fermi surfaces. Based on the three-band 
tight binding model, we obtain the specific heat, superfluid density, Knight shift and spin relaxation rate and 
find that all these properties at low temperature (T <C T c ) show powerlaw behaviors and are consistent avail¬ 
able experiments. Particularly, the superfluid density determined by the p z -wave pairing state in this quasi-one 
dimensional system is anisotropic: the in-plane superfluid density varies as Apy ~ T but the out-plane one 
varies as A p± ~ T 3 at low temperature. The anisotropic upper critical field reported in experiment is consis¬ 
tent with the S z = 0 (i.e., (fF + Ft)) Pz-wave pairing state. We also suggest the phase-sensitive dc-SQUID 
measurements to pin down the triplet p z -wave pairing state. 

PACS numbers: 74.20.-z, 74.25.-q, 74.70.-b, 74.20.Rp 


Recently, superconductivity in chromium-based materials 
has been revealed, attracting a lot of research interests 1-5 . 
For CrAs, the superconductivity is achieved by suppressing 
the magnetic order with pressure 1 - 2 . The close proximity 
of superconductivity to an helimagnetic order suggests an 
unconventional pairing mechanism 1 - 2 - 6 - 7 . The discovery of 
superconductivity in quasi-one dimensional(QlD) materials 
A2Cr 3 As 3 (A=K, Rb,Cs) 3 5 is especially interesting due to the 
rarity of Q 1 D superconductors. So far, the highest T c ~ 6.IK 
can been achieved in K2CT3AS3 at ambient pressure 3 . In nor¬ 
mal states, these materials possess large electronic specific- 
heat coefficients 3 and show non-Fermi liquid transport behav¬ 
ior. Particularly, K2CY3AS3 exhibits the exotic properties of 
Tomonaga-Luttinger liquids in NMR measurements 8 . In su¬ 
perconducting states, many unconventional superconducting 
properties have been experimentally observed, including lin¬ 
early temperature dependent penetration depth at T <C T c 9 , 
the absence of Hebel-Slichter coherence peak in I/T1 8 and the 
extremely large anisotropic upper critical field H c2 ( 0 ) 3-5 ’ 10-12 . 
All these suggest the existence of line nodes and possible spin 
triplet pairing in these Q 1 D superconductors. Furthermore, 
Raman scattering measurements suggest that electron phonon 
coupling is rather weak and magnetic fluctuations are coupled 
to the electronic structure via the lattice 13 . 

Theoretical calculations show that 3 d orbitals of Cr domi¬ 
nate the Fermi surfaces (FSs) consisting of two Q 1 D a and 
8 FSs and one three-dimensional 7 FS 14,15 . These Q 1 D 
materials possess strong frustrated magnetic fluctuations and 
are nearby a novel in-out co-planar magnetic ground state 15 . 
Moreover, the band structure near the Fermi level can be cap¬ 
tured by a minimum three-band tight binding model based 
on the A! x (d z 2) and E'(d xy ,d x 2_ y 2) molecular orbitals 16-18 . 
In a recent paper 17 , we adopted combined standard random 
phase approximation (RPA) approach to study the multi¬ 
orbital Hubbard-Hund model in the weak coupling limit and 
mean-field approach to study the t-J model in the strong cou¬ 
pling limit to investigate the pairing symmetry. Both ap¬ 


proaches consistently yield the triplet p,-wave pairing as the 
leading pairing symmetry for physically realistic parameters. 
The triplet pairing is driven by the ferromagnetic fluctuations 
within the sublattice. When considering spin-orbit coupling, 
the S z = 0 component (fF + Ft) slightly wins over the 
S z = ±1 ones (tt, 44 )- 

In this paper, we investigate the experimental observ¬ 
able properties of the triplet p~- wave pairing state in Q 1 D 
A2CT3AS3 materials, characterized by the line nodes on the 
k z = 0 plane on the Fermi surfaces. Based on the three-band 
tight binding model, we obtain the specific heat, superfluid 
density, Knight shift and spin relaxation rate and find that all 
these properties at low temperature show power-law behav¬ 
iors and are consistent available experiments. To be specific, 
when T <C T c , the specific heat varies as C v ~ T 2 , the 
Knight shift along the d-vector (|| z in this system) decays 
linearly but the other components remain almost unchanged 
and the spin relaxation rate varies as 1 /T\T ~ T 2 . Particu¬ 
larly, the superfluid density determined by the p z - wave pairing 
state in this Q 1 D system is anisotropic: the in-plane super¬ 
fluid density varies as py ~ T but the out-plane one varies as 
p± ~ T 3 at low temperature. The phase change along z di¬ 
rection in p z - wave state can be justified in the phase-sensitive 
measurements. The observed anisotropic upper critical field 
can be explained by the p 2 -wave pairing state with d || z (i.e. 

(tl + It))- 

Model: As the states near the Fermi level are contributed 
by three bands in K2Q3AS3, the minimum model to cap¬ 
ture the main physics is given by a three-band tight binding 
model 17 . Neglecting spin-orbit coupling (SOC), the tight- 
binding Hamiltonian is given by, 

J?TB = X^ t ( k ) /l ( k )^( k )> (D 

k 

where t/^k) = [cj t (k), 4 f (k), 4 t (k), c f 4 (k), 4 ; (k), 4 ; (k)]. 
The orbital index v = 1 , 2 , 3 represent the d z 2 for 1 , the d xy 
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for 2, and the d x 2_ y 2 for 3, respectively. The matrix h( k) has 
been given in Ref. 17. When considering the superconducting 
pair, the mean field Hamiltonian can be written as, 

H SC = ^^(k^sctk^k), 

k 

( h( k) A(k) \ 

%c( k ) — ^ A t(k) —/i*(—k) )’ (2) 

where ^(k) = [^(k), il> T (— k)]. Here we only consider the 
intraorbital pair, the pairing term A(k) is given by, 

A(k > = ( A«(k! £«(k) ) ■ (3 > 

/ A^'(k) \ 

A^'(k) = ;/ 2 (k) 1(4) 

V A^(k)/ 

In the following, unless otherwise specified, we shall con¬ 
sider the p 2 -wave spin triplet state with d = Asink z z, 
which means A^(k) = A^(k) = 0,A^(k) = A^(k) = 
A u sink z . Due to the existence of line nodes in the k z = 0 
plane, the quasi-particles exhibit linear density of states at low 
energies 17 . 

Diagonalizing the BdG matrix hsc( k) in Eq.(2), we obtain 
the eigenvalue E^ n as the band energy and the eigenvector 
</> n k as the eigen state. In the following calculations at fi¬ 
nite temperature T , the T —dependence of the gap amplitudes 
is assumed as A„(f) = A°<5(f). Here A° 2 _ j/2 = A° y = 
4A ° 2 = 0.83 meV are obtained in Ref.17, and S(t) is the 
normalized BCS gap at the reduced temperature t. = T/T c 19 . 
We have also performed calculations using the temperature 
dependent gap from solving the t-J model and found simi¬ 
lar temperature-dependence of the following experimental ob¬ 
servables at sufficiently low temperature T <C T c . 

Specific heat: The specific heat of this system is given by. 



FIG. 1: (color online) The electronic specific heat C v as a function 
of temperature for the p z -wave state. The inset shows the specific 
heat at low temperature. 

On the other hand, the temperature dependence of the super¬ 
fluid density obtained from the muon-spin relaxation mea¬ 
surements fits well to an isotropic s-wave character for the 
superconducting gap 21 . However, from our point of view, 
such conflicting experimental behaviors may be attributed to 
the anisotropic superfluid density of the Q1D p z -wave pairing 
state here. Actually, our results suggest that at low tempera¬ 
ture, while the in-plane superfluid density scales with T, the 
out-plane one scales with T 3 (which cannot be easily distin¬ 
guished from exponential function for T/T c « 1). 

The linear response of the system to an external magnetic 
field is detailed in the supplementary materials. In the super¬ 
conducting state, the superfluid density p is proportional to 
the response kernel I\ (q —> 0 , oj = 0 ), with p = p p + pd and 
K(q,uj) = K p (q,tu) + K d (q,co). In A 2 Cr 3 As 3 , the param¬ 
agnetic part p p and diamagnetic part pd read. 


C v 


1 

2N 


A’kn 


dn F (Ek n ) 

dT 


(5) 


where riF{Ek n ) = e E k „/l B T +1 is the Fermi-Dirac distribu¬ 
tion function. Here dnF ^ kri ' ) ca n be further evaluated as, 

dn F Ek n 1 dEk n 

Hr = "* (Ek “ )[1 - n *( J *“>](j£53 - kHfHtF '- <<s) 

The calculated specific heat C v is shown in Fig.l. Near T c , 
C v jumps, reflecting the superconducting phase transition. At 
low temperatures, C v oc T 2 (inset in Fig.l), reflecting the line 
gap nodes in the p,-wave state, consistent with experiments 20 . 

Superfluid density: The temperature-dependence of the su¬ 
perfluid density of the system takes on conflicting character¬ 
istics for different available experiments. On the one hand, 
through a tunnel diode oscillator, measurements of the tem¬ 
perature dependent penetration depth were done in K 2 Cr 3 As 3 
and a linear relationship was found at low temperature, which 
suggests the existence of line nodes in superconducting gap 0 . 


_ss 1 I jt pp/Ui 12 tt-F (77km) UF^Ekn) 

e, = W - £k „, _ Ek „ -• 

k mn 

P7 = ^El4k^(k)<MV(£k„). (7) 

k n 

Here s = x, y, z and Ff( k) and F^.(k) are given by, 

*?( k) = 

F? s ( k) = 


/s (k) 

° ^ 

(8) 

0 

-/r(- k) ) 5 

ft( k) 

0 ^ 

(9) 

0 

-f%(- k) ) ’ 


where ff (k) and f£ s ( k) are given in the supplementary ma¬ 
terials. For the intraband contribution, the Linhard function 
should be replaced with dn gj^ km ' > ■ 

The temperature dependent normalized superfluid density 
p is shown in Fig. 2. We find that p is anisotropic, which 
is caused by the p 2 -wave pairing state in this Q1D system. 
Firstly, the superfluid density along the z-axis is much larger 
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FIG. 2: (color online) The normalized superfluid density p ss (s = 
x, y, z) as a function of temperature for the p z -wave state. Due to 
the quasi-one dimensionality of the system, the absolute superfluid 
density p zz ( 0) : p xx ^ yv (0) « 20. The inset shows the superfluid 
density at low temperature. 


FIG. 3: (color online) The spin relaxation rate 1/TiT as a function 
of t£- for thep--wave state. Inset shows 1/T\T at low temperature. 


In general the spin susceptibility is defined as, 


than that in the xy -plane due to the Q1D characteristic of the 
systems. Secondly, while the superfluid density A p = p(T ) — 
p(0) in the xy -plane varies as A p» ~ T, the out-plane super¬ 
fluid density A p± varies slowly at low temperature, which 
seems at first glance like exponentially. However, further cal¬ 
culations show that the out-plane superfluid density varies as 
Ap± ~ T 3 , which cannot be easily distinguished from the 
exponential function at low T. To understand the qualitatively 
different behaviors of p || and p± at low temperatures, we have 
made a thorough investigation in the Supplementary Informa¬ 
tion. Here, we underdraw the physics. As the diamagnetic 
part pd is almost temperature independent, we can focus on 
the paramagnetic part p p , which represents for the consuming 
of the superfluid density by the nodal quasi-particle excita¬ 
tions aroused by the paramagnetic current operator J p . At low 
temperatures, these excitations mainly take place near k z = 0 
for the p z -wave pairing. For these small k z , the in-plane cur¬ 
rent is J v x j y ~ gi{k x , k y )( 1 — which leads to Apy ~ T. 
However, the out-plane current is ~ k z g 2 (k Xl k y ), which 
goes to zero for small k z and the aroused quasi-particle exci¬ 
tations are strongly suppressed, which leads to A p± ~ T 3 . 
Such low-temperature power-law behaviors of the superfluid 
density are the consequences of the polar states with an equa¬ 
torial line of nodes 22 2 '. The obtained anisotropic superfluid 
density for p.-wave state seems to be consistent with experi¬ 
ments, where both linear and seemingly exponential tempera¬ 
ture dependence of superfluid density were observed. 

NMR: Recently, the nuclear magnetic resonance(NMR) 
measurements suggest unconventional nature of superconduc¬ 
tivity, reflected in the absence of the Hebel-Slichter coherence 
peak below T c 8 . Moreover, the temperature dependence of 
1 j r l \ below T c shows power-law behavior. Now we investi¬ 
gate the spin-relaxation rate l/T), as well as the Knight shift 
K for the p 2 -wave triplet pairing state. 


rP 

= / (T T St(q,T)Sl(-q,0))e^ T dT,m 

Jo 

where .S'* is the t component of spin operator for orbital //. 
The NMR spin-relaxation rate reads, 
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AN 2 
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Y, A(k' -k)\<t>l k S s M 2 

,mns 

E=Ek rn &(Ekm T/i c / n ), (11) 


where the geometrical structure factor A(q) has been set to be 
1 in our calculation for simplicity. Here the spin matrices are 
S x , z = 773 (g> a x>z ® I 0 and S y = rj 0 ® cr y ( 8 ) I 0 , where the 
Pauli matrices a and 77 act in the spin and Nambu space and 
the 3x3 identity Io acts in the orbital space. 

Fig. 3 shows the spin relaxation rate \/T{T for the p 2 -wave 
state. The most distinctive feature of 1 /T\T is that there is no 
Hebel-Slichter coherence peak below T c . At low temperature, 
1/TiT shows powerlaw behavior, i.e., 1/T\ ~ T 3 , because 
of the line nodes in superconducting gap. Such results are 
consistent with the experimental data although the power law 
exponent is slightly different from the experimental one 8 . 

The Knight shift reads. 


Kss cx 5>-(0,0) 


1 

AN 


T, 

k mn 


TlF^Ekm) 


^*km 


riF(Ekn) 

E\tn 


( 12 ) 


The Knight shift experiment can be used to distinguish be¬ 
tween spin-singlet and spin-triplet pairings. Further more, for 
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FIG. 4: (color online) The Knight shift K S s(s = x, y, z) as a func¬ 
tion of temperature for the p z -wave state, (a) d = Asink,z and the 
corresponding parameter is A 44 = A 44 , (b) d = A smfc z x and the 
corresponding parameter is A 44 = —A”. 


the triplet pairing, it can identify the multi-component order 
parameter through the d-vector defined as, 


d k = ^ 


Iktt 


,A kti ]. (13) 



(b) (c) 



The obtained Knight shifts K ss for the d = Asink z z and 
d = A sink z x pairing states are shown in Fig.4(a) and (b), 
respectively. In the superconducting state, the component of 
Knight shift along d is strongly suppressed while the other 
components almost unchange with the decreasing of temper¬ 
ature. For the d = Asink z z state obtained in Ref. 17, the 
spin orientation is confined in the xy- plane. Therefore, the 
in-plane magnetic field can polarize the spin without breaking 
the Cooper pair. On the contrary, the out-plane magnetic field 
need to break the Cooper pairing in order to cause spin split¬ 
ting. Therefore, this excitation is linearly suppressed at low 
temperature due to the existence of line nodes in supercon¬ 
ducting state. Similarly, we can understand the Knight shift 
for the d = A sink z x state. 

Phase sensitive measurements: The triplet p z - wave pair¬ 
ing predicted here can be detected by the dc SQUID, a 
phase-sensitive device which has been adopted in determin¬ 
ing the pairing symmetries of such superconducting sys¬ 
tems as the cuprates 24-26 , the Sr 2 Ru 04 27 and others 2 *. The 
proposed configuration is shown in Fig. 5(a), where two 
superconductor-normal metal-superconductor (SNS) Joseph- 
son tunneling junctions are formed on the two opposite edges 
in the z direction of K 2 Cr 3 As 3 , which are connected by a loop 
of a conventional s-wave superconductor, forming a bimetal¬ 
lic ring with a magnetic flux $ threading through the loop. 

As a result of the interference between the two branches of 
Josephson supercurrent, the maximum total supercurrent (the 
critical current) I c in the circuit modulates with >I> according 
to 


Ic($) = 2/ 0 


cos 


$ . </>0 


(14) 


FIG. 5: (color online) Experimental geometry for a SQUID phase 
sensitive probe and interference patterns of a SQUID, (b) s-wave 
case. (c)p z -wave case. 


Here Iq is the critical current of one Josephson junction, <f>o = 
h/2e is the basic flux quantum, 4>o is the phase shift along z 
direction of K 2 C 13 AS 3 . 

If K 2 CT 3 AS 3 is a conventional superconductor, there is no 
phase change along z direction so that <fi 0 = 0 and the maxi¬ 
mum current should show a maximum for > 1 ) = n<l>o, as shown 
in Fig. 5(b). However, for p z - wave state there a n phase shift 
along z direction, that is, Qq = 7 r. The corresponding current 
should show a maximum for $ = (n + 4 )<l>o- as shown in 
Fig. 5(c). The phase change obtained in this experiment can 
provide a verdict for the p,-wave pairing state. 

Upper critical magnetic field: One of the most striking 
features of A 2 CT 3 AS 3 superconductors is that the upper crit¬ 
ical field H c2 is very high and severely exceeds the Pauli 
limit 3-510-12 . Furthermore, the H c2 along the chains shows 
paramagnetic-limited behavior but the H c 2 perpendicular to 
the chains does not, which results in a smaller H c2 along z 
axis at 0 K. The observed high H c2 can hardly be attributed 
to the spin orbital effect but may suggest triplet pairing. Ac¬ 
tually, the high anisotropic H c2 observed in these materials is 
consistent with the proposed p~-wave spin triplet state with 
d = Asink z z (i.e. (fj, + j,j')). The magnetic field par¬ 
allel to d(z direction) will cause Zeeman splitting and break 
the Cooper pair, which leads to the paramagnetic-limited be¬ 
havior. However, the in-plane magnetic field will not cause 
pair breaking because this pairing can be equally looked upon 
as the iteration of two equal spin pairings with opposite spin 
directions in the x, y-plane. Thus, there is no paramagnetic 
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suppression of superconductivity for in-plane magnetic field. 
The extraordinary crossover of H C 2 (T) curves in K^Cr-j As :j is 
rather similar to that of a heavy fermion system UPt 3 29 . The 
crossover in the latter can hardly be explained by a spin sin¬ 
glet state but can be well explained by a spin triplet state with 
SOC 30,31 . Therefore, the crossover in K 2 Cr 3 As 3 may be simi¬ 
larly understood to be caused by the p z -wave spin triplet state 
with SOC. Further detailed calculations are needed. 

In conclusion, we have investigated the experimental con¬ 
sequences of the p z -wave spin triplet superconductivity in 
A 2 Cr 3 As 3 . Based on the three-band tight binding model, we 
obtain the specific heat, superfluid density, Knight shift and 
spin relaxation rate and find that all these properties at low 
temperature show powlaw temperature-dependences and are 
consistent with available experiments. Due to the existence of 
line nodes in k z = 0 plane, the quasi-particles exhibit linear 
density of states at low energies. The specific heat exhibits 
a quadratic temperature dependence at T <C T c . The p z - 
wave pairing state leads to an anisotropic superfluid density 
in the Q1D materials: the in-plane superfluid density varies 


as Apy ~ T but the out-plane one varies as Ap^ ~ T 3 at 
low temperature, which is consistent with experimental data. 
In the spin relaxation rate, there is no Hebel-Slichter coher¬ 
ence peak below T c and 1 /1\ varies as T 3 at low temperature. 
In our obtained superconducting state with d = Asink z z 
(which means a (f| + it) pairing), the out-plane Knight 
shift is linearly suppressed at low temperature but the in¬ 
plane one remains almost unchanged. The phase stmcture in 
the p ; .-wave state can be justified in the phase-sensitive dc- 
SQUID measurements. Finally, the experimentally observed 
anisotropic upper critical field can be explained by the p z - 
wave pairing state with d j| z obtained here. 
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Supplementary materials for “Experimental consequences of p z -wave spin triplet superconductivity 

in A 2 Cr 3 As 3 ” 


1. The linear response of the system to an external magnetic field 

In the presence of a weak external magnetic field taken as a perturbation, the coupling of the system to this field can be 
introduced via the Peierls substitution c\ arT c :j p a —> c\ aa e ^ Adl Cjp a in the tight-binding Hamiltonian (1), where A is the vector 
potential. This Hamiltonian can be expanded up to the second-order terms with respect to A as 

H( A) = H( 0) - ^[J p (f)A(i) - l -K d (i) A 2 (i)], (SI) 

where J p is the paramagnetic current. Then the current J s (i) (s = x, y, z ) is given by. 

Mi) = -jrrl = (S2) 

dA s (i) 

In momentum space, the total current can be further written as a response of the system to the vector potential, 

Mq) = + K i(q)]Mq)- (S3) 

t 

Here K p st and K d t are the paramagnetic and diamagnetic response kernels, respectively. The diamagnetic response kernel in 
this system reads. 


K d s M) = £>t(k)/* (k + q)t/;(k + q)>, (S4) 

k 

with f d t (k + q) = |k=k+q- The paramagnetic response kernel is obtained through the current-current correlation function 

in linear response, 

/ OO 

_ t'){[J p {q, t), J? (- q, t')])dt, (S5) 

-OO 

where J p (q) = '0^(k)/|’(k + q)V’(k + q) and /f(k + q) = d gl k) [k=k+q- The corresponding Matsubara Green function 

is, 


K st(QA^n)= [ (TVT p (q,T)J t p (-g,0))e l ™'‘dr (S6) 

Jo 


with uj n — 1 . The retarded Green function can be obtained upon analytic continuation, K p t (q : u>) = K p t {q , iuj n —> oj + id). 


2. The current operators 


The paramagnetic current operator is given by J p (q) = ]T k -i/A (k)/ p (k + q)^(k + q) and the diamagnetic response kernel is 
K d sM) = Ek(^(k)/ S d t (k + q)V(k + q)). In the following, we provide the explicit formula of the above / p and f d t matrices. 
Let x = ^k x ao , y = \k y aa and z = \k z c$, the matrix elements / p (k) and f d z ( k) are. 


fP 
J zll 

fP 
Jz 22 
fP 

Jz33 


f d 

J zz 11 
f d 

J zz22 
f d 

J zz 33 


—2cos^sin22(2cos2j/ + 4coscccos?/) — 2cas 1 3 1 sin2z — 4cos].\sin4z — 6coslgSin6z, 
—2coS 2 J±sin2z — 4coS?4sm4;r 

cp fP _ fP _ fP _ n 
J z22 J J 312 ~ Jz 13 — J z23 ~ u > 

— 2c%s y lcos2z(2cos2y + 4cosxcosy ) — 2cqS*2Cos2z — ?>6qsV^cos4z — 18cQS*gCos6z, 
—2cqS^2Cos2^ — 8cqS 22 cos4z 

fd fd _ fd _ fd _ 

J zz22i Jzzl2 Jzz 13 Jzz23 ^ 


(S7) 


(S8) 
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The matrix elements f?( k) and f yy ( k) are, 

fy U = -2a 0 (sll + 2s 1 y 1 z cos2z){sva2y + cosxsim/) 

fyi 2 = 2iaosJ 2 cos2y + iaos y 2 cosycosx + V3aoS2yCOsysmx 

f yl3 = —2aoS2ySin2y + i\Z3aoSi y sim/sinx + aos^sinycosx 

/y22 = 2cio s iij / siia2z/ — — (to( s iiy T 3s22y) cos 2 :s i n 2/ 

fy23 = ^Y a o(s'n y - sii !/ )sinxcost/ + 2ia 0 sllyCOs2y - 2ia 0 sH y cosxcosy 

fy 33 = -2a 0 sfi y sm2y - ^a 0 (3s^ y + s^ y )cosxsiny, 


fyyil = -2al(s l y \ + 2s l y Z cos2z)(cos2y + -cosxcosy) 

fyyl2 = 

fyy 13 = 

fyy22 = 


— 2ials y 2 sm2y — - ia^s y 2 sinycosx — -^-OgS^sini/sina; 
— 2ciQS2yCOs2y + i^-aQs\ 2 cosysinx + -ags^cosj/cosa; 
-2ogs?; w cos2y - + Ss 2 ^)cosxcosy 


/; 


J/l/23 


rd 

Jyy33 


2 0 S 11y c 

V3 

4 


—— ao( s iiy ~ s 22i/) s i na ' s i n y — 2?'aQS 22 y sin2y + iaos 2 2 y cosa;sir] 


—2agS 22 y c °s2j/ - ^a^(3s^ y + si|y>cosxcosy 


The matrix elements /J?(k) and / xx (k) are, 

/ xll = -2v / 3a 0 (s y 1 + 2s**cos2z)sinxcos2/ 

/ xl2 = —%/3iaosJ 2 sinysinx + 3aoS2 2 sim/cosx 
/ xl3 = — 3?'aos} 2 cos'(/cosx + v^aos^yCosysi 11 ^ 

fx22 = ~~^ a o{ s lly + 3S22y)sinXCOSt/ 

/x23 = ^«o(s?i y - s 22i/)cosxsin7/ + 2v / 3*a 0 s 2 | y sinxsinj/ 

/x33 = -^ a o(3s?r y + s|y sinxcosy, 


f d 

J xxll 
fd 

J xxl2 

fd 

J xxl3 

fd 

J xx22 
fd. 

J xx23 

fd 

J xx33 


—Saj^s* 1 + 2 s y ^cos 22 :)cosxcosy 


3 . 2 12 . • 3 %/3 


2 „ 12 r . 


— -za 0 Si y smycosx- —a 0 s 2y smysmxx 

3 y/3 .212 • , 3 2 12 

ta 0 5 ly cosj/smx + -a 0 s 2y cosycosx 

-| a o( s lly + 3s| 2 y) cosxcosy 

--^- a o( s iiy - slly)smxsiny + 3 ialsll y cosxsmy 
-^ a o( 3s ii y + sll y )cosxcosy 


(S9) 


(S10) 


(Sll) 


(S12) 
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3. Superfluid density at low temperature for p z -wave pairing state 


Here we investigate the temperature-dependence of the superfluid density of the p z -wave pairing state in this Q1D system at 
low temperature T CT c . The diamagnetic part of the superfluid density is almost temperature independent at low temperature, 
and therefore we shall focus on the paramagnetic part, which represents for the consuming of the superfluid density through the 
nodal quasi-particle excitation near k z = 0. At low temperature, only the intra-band quasi-particle excitation is important, which 
dictates us to simplify the paramagnetic part in Eq.(7) of the main text as, 



(S13) 



(S14) 


km 


At low temperature T <C T c « 0( A„), the derivative term = — l/[4fcsT cosh 2 (E^ m /2ksT)] is obviously non¬ 

zero only for those momenta km with the band energy £k m satisfying > |i?k m | ~ |Akm| oc \k z \, which occupy the 
region near the line gap nodes at k z = 0. Thus the momentum summation in Eq.(S13) and Eq.(S14) is dominantly contributed 
from the region near k z = 0. At that region, from Eq.(8) in the main text and Eq.(S7), Eq.(S9) and Eq.(Sll), we can verify 
that ^Lk^/ S ( k )^k ~ gi{k x ,k y ){a - k 2 z ) and ^ k ^I( k )4k ~ g 2 {k x ,k y )k z . With these approximations, the in-plane 
superfluid density pp X ^ vv can be further estimated as. 



(S15) 


The out-plane superfluid density p^ z is. 



fS 16) 


Note that on the above, we have used the relation N(E) ex \E\ for the low energy excitations near the line gap nodes. 











